Hypoxia treatment and resistance training alters microRNA profiling in rats skeletal muscle.

MicroRNAs (miRNAs) may play a crucial regulatory role in the process of muscle atrophy induced by high-altitude hypoxia and its amelioration through resistance training. However, research in this aspect is still lacking. Therefore, this study aimed to employ miRNA microarray analysis to investigate the expression profile of miRNAs in skeletal muscle from an animal model of hypoxia-induced muscle atrophy and resistance training aimed at mitigating muscle atrophy. The study utilized a simulated hypoxic environment (oxygen concentration at 11.2%) to induce muscle atrophy and established a rat model of resistance training using ladder climbing, with a total intervention period of 4 weeks. The miRNA expression profile revealed 9 differentially expressed miRNAs influenced by hypoxia (e.g., miR-341, miR-32-5p, miR-465-5p) and 14 differentially expressed miRNAs influenced by resistance training under hypoxic conditions (e.g., miR-338-5p, miR-203a-3p, miR-92b-3p) (∣log2(FC)∣ ≥ 1.5, p < 0.05). The differentially expressed miRNAs were found to target genes involved in muscle protein synthesis and degradation (such as Utrn, mdm2, eIF4E), biological processes (such as negative regulation of transcription from RNA polymerase II promoter, regulation of transcription, DNA-dependent), and signaling pathways (such as Wnt signaling pathway, MAPK signaling pathway, ubiquitin-mediated proteolysis, mTOR signaling pathway). This study provides a foundation for understanding and further exploring the molecular mechanisms underlying hypoxia-induced rats muscle atrophy and the mitigation of atrophy through resistance training.

Genome-wide identification and expression analysis of ARF gene family in embryonic development of Korean pine (Pinus koraiensis).

BACKGROUND: The Auxin Responsive Factor (ARF) family plays a crucial role in mediating auxin signal transduction and is vital for plant growth and development. However, the function of ARF genes in Korean pine (Pinus koraiensis), a conifer species of significant economic value, remains unclear. RESULTS: This study utilized the whole genome of Korean pine to conduct bioinformatics analysis, resulting in the identification of 13 ARF genes. A phylogenetic analysis revealed that these 13 PkorARF genes can be classified into 4 subfamilies, indicating the presence of conserved structural characteristics within each subfamily. Protein interaction prediction indicated that Pkor01G00962.1 and Pkor07G00704.1 may have a significant role in regulating plant growth and development as core components of the PkorARFs family. Additionally, the analysis of RNA-seq and RT-qPCR expression patterns suggested that PkorARF genes play a crucial role in the development process of Korean pine. CONCLUSION: Pkor01G00962.1 and Pkor07G00704.1, which are core genes of the PkorARFs family, play a potentially crucial role in regulating the fertilization and developmental process of Korean pine. This study provides a valuable reference for investigating the molecular mechanism of embryonic development in Korean pine and establishes a foundation for cultivating high-quality Korean pine.

Inhibition of signaling protein ERN1 increases the sensitivity of serine synthesis gene expressions to glucose and glutamine deprivations in U87MG glioblastoma cells.

Objective. Glucose and glutamine supply as well as serine synthesis and endoplasmic reticulum (ER) stress are important factors of glioblastoma growth. Previous studies showed that the knockdown of ERN1 (ER to nucleus signaling 1) suppressed glioblastoma cell proliferation and modified the sensitivity of numerous gene expressions to nutrient deprivations. The present study is aimed to investigate the impact of glucose and glutamine deprivations on the expression of serine synthesis genes in U87MG glioblastoma cells in relation to ERN1 knockdown with the intent to reveal the role of ERN1 signaling pathway on the ER stress-dependent regulation of these gene expressions. Clarification of the regulatory mechanisms of serine synthesis is a great significance for glioblastoma therapy. Methods. The control U87MG glioblastoma cells (transfected by empty vector) and ERN1 knockdown cells (transfected by dominant-negative ERN1) were exposed under glucose and glutamine deprivation conditions for 16 h. RNA was extracted from cells and reverse transcribed. The expression level of PHGDH (phosphoglycerate dehydrogenase), PSAT1 (phosphoserine amino-transferase 1), PSPH (phosphoserine phosphatase), ATF4 (activating transcription factor 4), and SHMT1 (serine hydroxymethyltransferase 1) genes was studied by real-time qPCR and normalized to ACTB. Results. It was found that the expression level of genes responsible for serine synthesis such as PHGDH, PSAT1, PSPH, and transcription factor ATF4 was up-regulated in U87MG glioblastoma cells under glucose and glutamine deprivations. Furthermore, inhibition of ERN1 significantly enhances the impact of glucose and especially glutamine deprivations on these gene expressions. At the same time, the expression of the SHMT1 gene, which is responsible for serine conversion to glycine, was down-regulated in both nutrient deprivation conditions with more significant changes in ERN1 knockdown glioblastoma cells. Conclusion. Taken together, the results of present study indicate that the expression of genes responsible for serine synthesis is sensitive to glucose and glutamine deprivations in gene-specific manner and that suppression of ERN1 signaling significantly modifies the impact of both glucose and glutamine deprivations on PHGDH, PSAT1, PSPH, ATF4, and SHMT1 gene expressions and reflects the ERN1-mediated genome reprograming introduced by nutrient deprivation condition.

Skin in the game: a review of single-cell and spatial transcriptomics in dermatological research.

Single-cell RNA sequencing (scRNA-seq) and spatial transcriptomics (ST) are two emerging research technologies that uniquely characterize gene expression microenvironments on a cellular or subcellular level. The skin, a clinically accessible tissue composed of diverse, essential cell populations, serves as an ideal target for these high-resolution investigative approaches. Using these tools, researchers are assembling a compendium of data and discoveries in healthy skin as well as a range of dermatologic pathophysiologies, including atopic dermatitis, psoriasis, and cutaneous malignancies. The ongoing advancement of single-cell approaches, coupled with anticipated decreases in cost with increased adoption, will reshape dermatologic research, profoundly influencing disease characterization, prognosis, and ultimately clinical practice.

Association of gene polymorphisms and the decreased expression of long non-coding RNA LOC553103 with rheumatoid arthritis.

BACKGROUND: Long non-coding RNAs (lncRNAs) are involved in many key bioprocesses, including the occurrence and development of rheumatoid arthritis (RA). We aimed to analyze the association of genetic variants of long non-coding RNA LOC553103 and its peripheral blood mononuclear cells (PBMC) expression with RA. METHODS: We enrolled 457 RA patients and 551 healthy controls and conducted a case-control study to analyze the relationship between LOC553103 gene rs272879 and the susceptibility of RA by TaqMan single nucleotide polymorphism genotyping. Among them, we sampled 92 cases and 92 controls, respectively, to detect the PBMC level of LOC553103 using quantitative real-time polymerase chain reaction technology. We explored the association between LOC553103 rs272879 and its PBMC expression levels in 71 RA patients. Mann-Whitney, Chi-square, and Spearman correlation analysis were used for statistical analysis and P-value <.05 was considered statistically significant. RESULTS: The genotype frequency of LOC553103 rs272879 CC was increased, and CG was decreased in RA patients compared to the control group (χ2 = 6.772, P = .034). The LOC553103 expression level in PBMC of RA patients was downregulated compared to healthy control (Z = -4.497, P < .001). Moreover, negative correlations were observed between the PBMC level of LOC553103 and erythrocyte sedimentation rate (rs = -0.262, P = .018), white blood cell count (rs = -0.382, P = .004), platelet (rs = -0.293, P = .030), and disease activity score in 28 joints (rs = -0.271, P = .016) in RA patients. CONCLUSIONS: This study provides the first evidence supporting an association between LOC553103 gene polymorphisms and susceptibility of RA and a relationship of PBMC level of LOC553103 with clinical manifestations and laboratory indicators of RA patients.

A General Protocol for Gene Expression Analysis in Chemically Mutant Coffee Plants in Response to Rust (Hemileia vastatrix Berk & Broome) Infection.

Coffea spp. is the source of one of the most widely consumed beverages in the world. However, the cultivation of this crop is threatened by Hemileia vastatrix Berk & Broome, a fungal disease, which reduces the productivity and can cause significant economic losses. In this protocol, coffee leaf segment derived from a chemical mutagenesis process are inoculated with uredospores of the pathogen. Subsequently, the gene expression changes are analyzed over the time (0, 5, 24, 48, and 120 h) using quantitative real-time polymerase chain reaction (RT-qPCR). The procedures and example data are presented for expression analysis in the CaWRKY1 gene. This procedure can be applied for quantitative analysis of other genes of interest to coffee breeders and scientists for elucidating the molecular mechanisms involved in the interaction between the plant and pathogen, potentially leading to the development of more efficient approaches for managing this disease.

Transcriptomic Data Analysis Using the Galaxy Platform: Coffee (Coffea arabica L.) Flowers as Example.

Coffee, an important agricultural product for tropical producing countries, is facing challenges due to climate change, including periods of drought, irregular rain distribution, and high temperatures. These changes result in plant water stress, leading to significant losses in coffee productivity and quality. Understanding the processes that affect coffee flowering is crucial for improving productivity and quality. In this chapter, we describe a protocol for transcriptome analysis using available Internet software, mainly in the Galaxy Platform, using RNA-Seq data from flowers collected from different parts of the coffee tree. The methods presented in this chapter provide a comprehensive protocol for transcriptome analysis of differentially expressed genes from flowers of coffee plant. This knowledge can be utilized in coffee genetic improvement programs, particularly in the selection of cultivars that are tolerant to water deficit.

Quantitative Estimation of Promoter Activity in Cannabis sativa Using Agroinfiltration-Based Transient Gene Expression.

To properly assess promoter activity, which is critical for understanding biosynthetic pathways in different plant species, we use agroinfiltration-based transient gene expression assay. We compare the activity of several known promoters in Nicotiana benthamiana with their activity in Cannabis sativa (both hemp and medicinal cannabis), which has attracted much attention in recent years for its industrial, medicinal, and recreational properties. Here we describe an optimized protocol for transient expression in Cannabis combined with a ratiometric GUS reporter system that allows more accurate evaluation of promoter activity and reduces the effects of variable infiltration efficiency.

Agroinfiltration for Enhanced Transgene Expression in Coffee Leaves (Coffea arabica L.).

The Coffea spp. plant is a significant crop in Latin America, Africa, and Asia, and recent advances in genomics and transcriptomics have opened possibilities for studying candidate genes and introducing new desirable traits through genetic engineering. While stable transformation of coffee plants has been reported using various techniques, it is a time-consuming and laborious process. To overcome this, transient transformation methods have been developed, which avoid the limitations of stable transformation. This chapter describes an ex vitro protocol for transient expression using A. tumefaciens-mediated infiltration of coffee leaves, which could be used to produce coffee plants expressing desirable traits against biotic and abiotic stresses, genes controlling biochemical and physiological traits, as well as for gene editing through CRISPR/Cas9.

Transcriptional elongation control of hypoxic response.

The release of paused RNA polymerase II (RNAPII) from promoter-proximal regions is tightly controlled to ensure proper regulation of gene expression. The elongation factor PTEF-b is known to release paused RNAPII via phosphorylation of the RNAPII C-terminal domain by its cyclin-dependent kinase component, CDK9. However, the signal and stress-specific roles of the various RNAPII-associated macromolecular complexes containing PTEF-b/CDK9 are not yet clear. Here, we identify and characterize the CDK9 complex required for transcriptional response to hypoxia. Contrary to previous reports, our data indicate that a CDK9 complex containing BRD4 but not AFF1/4 is essential for this hypoxic stress response. We demonstrate that BRD4 bromodomains (BET) are dispensable for the release of paused RNAPII at hypoxia-activated genes and that BET inhibition by JQ1 is insufficient to impair hypoxic gene response. Mechanistically, we demonstrate that the C-terminal region of BRD4 is required for Polymerase-Associated Factor-1 Complex (PAF1C) recruitment to establish an elongation-competent RNAPII complex at hypoxia-responsive genes. PAF1C disruption using a small-molecule inhibitor (iPAF1C) impairs hypoxia-induced, BRD4-mediated RNAPII release. Together, our results provide insight into potentially targetable mechanisms that control the hypoxia-responsive transcriptional elongation.

Gene expression profiling in elderly patients with familial hypercholesterolemia with and without coronary heart disease.

BACKGROUND AND AIMS: Elderly familial hypercholesterolemia (FH) patients are at high risk of coronary heart disease (CHD) due to high cholesterol burden and late onset of effective cholesterol-lowering therapies. A subset of these individuals remains free from any CHD event, indicating the potential presence of protective factors. Identifying possible cardioprotective gene expression profiles could contribute to our understanding of CHD prevention and future preventive treatment. Therefore, this study aimed to investigate gene expression profiles in elderly event-free FH patients. METHODS: Expression of 773 genes was analysed using the Nanostring Metabolic Pathways Panel, in peripheral blood mononuclear cells (PBMCs) from FH patients ≥65 years without CHD (FH event-free, n = 44) and with CHD (FH CHD, n = 39), and from healthy controls ≥70 years (n = 39). RESULTS: None of the genes were differentially expressed between FH patients with and without CHD after adjusting for multiple testing. However, at nominal p < 0.05, we found 36 (5%) differentially expressed genes (DEGs) between the two FH groups, mainly related to lipid metabolism (e.g. higher expression of ABCA1 and ABCG1 in FH event-free) and immune responses (e.g. lower expression of STAT1 and STAT3 in FH event-free). When comparing FH patients to controls, the event-free group had fewer DEGs than the CHD group; 147 (19%) and 219 (28%) DEGs, respectively. CONCLUSIONS: Elderly event-free FH patients displayed a different PBMC gene expression profile compared to FH patients with CHD. Differences in gene expression compared to healthy controls were more pronounced in the CHD group, indicating a less atherogenic gene expression profile in event-free individuals. Overall, identification of cardioprotective factors could lead to future therapeutic targets.

Cis-regulatory control of transcriptional timing and noise in response to estrogen.

Cis-regulatory elements control transcription levels, temporal dynamics, and cell-cell variation or transcriptional noise. However, the combination of regulatory features that control these different attributes is not fully understood. Here, we used single-cell RNA-seq during an estrogen treatment time course and machine learning to identify predictors of expression timing and noise. We found that genes with multiple active enhancers exhibit faster temporal responses. We verified this finding by showing that manipulation of enhancer activity changes the temporal response of estrogen target genes. Analysis of transcriptional noise uncovered a relationship between promoter and enhancer activity, with active promoters associated with low noise and active enhancers linked to high noise. Finally, we observed that co-expression across single cells is an emergent property associated with chromatin looping, timing, and noise. Overall, our results indicate a fundamental tradeoff between a gene's ability to quickly respond to incoming signals and maintain low variation across cells.

Tissue-specific gene expression of genome-wide significant loci associated with major depressive disorder subtypes.

Major depressive disorder (MDD) is a clinically and genetically heterogeneous disorder. To reduce heterogeneity, large-scale genome-wide association studies have recently identified genome-wide significant loci associated with seven MDD subtypes. However, it was unclear in which tissues the genes near those loci are specifically expressed. We investigated whether genes related to specific MDD subtypes would be preferably expressed in a specific tissue. At 14 novel subtype-specific loci related to seven MDD subtypes-(1) non-atypical-like features MDD, (2) early-onset MDD, (3) recurrent MDD, (4) MDD with suicidal thoughts, (5) MDD without suicidal thoughts, (6) MDD with moderate impairment, and (7) postpartum depression, we investigated whether 22 genome-wide significant genetic variant-mapped genes were tissue-specifically expressed in brain, female reproductive, male specific, cardiovascular, gastrointestinal, or urinary tissues in the Genotype-Tissue Expression (GTEx) subjects (n ≤ 948). To confirm the tissue-specific expression in the GTEx, we used independent Human Protein Atlas (HPA) RNA-seq subjects (n ≤ 95). Of 22 genes, nine and five genes were tissue-specifically expressed in brain and female reproductive tissues, respectively (p < 2.27 × 10-3). RTN1, ERBB4, and AMIGO1 related to early-onset MDD, recurrent MDD, or MDD with suicidal thoughts were highly expressed in brain tissues (d = 1.19-2.71), while OAS1, LRRC9, DHRS7, PSMA5, SYPL2, and GULP1 related to non-atypical-like features MDD, early-onset MDD, MDD with suicidal thoughts, or postpartum depression were expressed at low levels in brain tissues (d = -0.17--1.48). DFNA5, CTBP2, PCNX4, SDCCAG8, and GULP1, which are related to early-onset MDD, MDD with moderate impairment, or postpartum depression, were highly expressed in female reproductive tissues (d = 0.80-2.08). Brain and female reproductive tissue-specific expression was confirmed in the HPA RNA-seq subjects. Our findings suggest that brain and female reproductive tissue-specific expression might contribute to the pathogenesis of MDD subtypes.

Early life adversity is associated with differential gene expression in immune cells: A cluster-based analysis across an acute psychosocial stressor.

Elucidating mechanisms by which early-life adversity (ELA) contributes to increased disease risk is important for mitigating adverse health outcomes. Prior work has found differences in immune cell gene expression related to inflammation and mitochondrial activity. Using a within-person between-group experimental design, we investigated differences in gene expression clusters across acute psychosocial stress and no-stress conditions. Participants were young adults (N = 29, aged 18 - 25 years, 62 % female, 47 % with a history of ELA). Gene expression was assessed in peripheral blood mononuclear cells collected at 8 blood draws spanning two 5-hour sessions (stress vs. no-stress) separated by a week, 4 across each session (number of observations = 221). We applied two unsupervised gene clustering methods - latent profile analysis (LPA) and weighted gene co-expression analysis (WGCNA) - to cluster genes with similar expression patterns across participants. LPA identified 11 clusters, 7 of which were significantly associated with ELA-status. WGCNA identified 5 clusters, 3 of which were significantly associated with ELA-status. LPA- and WGCNA-identified clusters were correlated, and all clusters were highly preserved across sessions and time. There was no significant effect of acute stress on cluster gene expression, but there was a significant effect of time, and significant differences by ELA-status. ELA-associated clusters related to RNA splicing/processing, inflammation, leukocyte differentiation and division, and mitochondrial activity were differentially expressed across time: ELA-exposed individuals showed decreased expression of these clusters at 90-minutes while controls showed increased expression. Our findings replicate previous work in this area and highlight additional mechanisms by which ELA may contribute to disease risk.

Identification of genes related to ribosomal proteins in colorectal cancer: exploring their potential as biomarkers, prognostic indicators, and therapeutic targets.

BACKGROUND: Colorectal cancer (CRC) ranks as the third most commonly diagnosed cancer in both females and males, underscoring the need for the identification of effective biomarkers. METHODS AND RESULTS: We assessed the expression levels of ribosomal proteins (RPs) at both mRNA and protein levels. Subsequently, leveraging the STRING database, we constructed a protein-protein interaction network and identified hub genes. The co-expression network of differentially expressed genes associated with CRC and their target hub RPs was constructed using the weighted gene co-expression network analysis algorithm. Gene ontology and molecular signatures database were conducted to gain insights into the biological roles of genes associated with the identified module. To confirm the results, the expression level of the candidate genes in the CRC samples compared to the adjacent healthy was evaluated by the RT-qPCR method. Our findings indicated that the genes related to RPs were predominantly enriched in biological processes associated with Myc Targets, Oxidative Phosphorylation, and cell proliferation. Also, results demonstrated that elevated levels of GRWD1, MCM5, IMP4, and RABEPK that related to RPs were associated with poor prognostic outcomes for CRC patients. Notably, IMP4 and RABEPK exhibited higher diagnostic value. Moreover, the expression of IMP4 and RABEPK showed a significant association with drug resistance using cancer cell line encyclopedia and genomics of drug sensitivity in cancer databases. Also, the results showed that the expression level of IMP4 and RABEPK in cancerous samples was significantly higher compared to the adjacent healthy ones. CONCLUSION: The general results of this study have shown that many genes related to RPs are increased in cancer and could be associated with the death rate of patients. We also highlighted the therapeutic and prognostic potentials of RPs genes in CRC.

Joint Bayesian estimation of cell dependence and gene associations in spatially resolved transcriptomic data.

Recent technologies such as spatial transcriptomics, enable the measurement of gene expressions at the single-cell level along with the spatial locations of these cells in the tissue. Spatial clustering of the cells provides valuable insights into the understanding of the functional organization of the tissue. However, most such clustering methods involve some dimension reduction that leads to a loss of the inherent dependency structure among genes at any spatial location in the tissue. This destroys valuable insights of gene co-expression patterns apart from possibly impacting spatial clustering performance. In spatial transcriptomics, the matrix-variate gene expression data, along with spatial coordinates of the single cells, provides information on both gene expression dependencies and cell spatial dependencies through its row and column covariances. In this work, we propose a joint Bayesian approach to simultaneously estimate these gene and spatial cell correlations. These estimates provide data summaries for downstream analyses. We illustrate our method with simulations and analysis of several real spatial transcriptomic datasets. Our work elucidates gene co-expression networks as well as clear spatial clustering patterns of the cells. Furthermore, our analysis reveals that downstream spatial-differential analysis may aid in the discovery of unknown cell types from known marker genes.

Gene expression analysis of drought tolerance and cuticular wax biosynthesis in diploid and tetraploid induced wallflowers.

Whole-genome doubling leads to cell reprogramming, upregulation of stress genes, and establishment of new pathways of drought stress responses in plants. This study investigated the molecular mechanisms of drought tolerance and cuticular wax characteristics in diploid and tetraploid-induced Erysimum cheiri. According to real-time PCR analysis, tetraploid induced wallflowers exhibited increased expression of several genes encoding transcription factors (TFs), including AREB1 and AREB3; the stress response genes RD29A and ERD1 under drought stress conditions. Furthermore, two cuticular wax biosynthetic pathway genes, CER1 and SHN1, were upregulated in tetraploid plants under drought conditions. Leaf morphological studies revealed that tetraploid leaves were covered with unique cuticular wax crystalloids, which produced a white fluffy appearance, while the diploid leaves were green and smooth. The greater content of epicuticular wax in tetraploid leaves than in diploid leaves can explain the decrease in cuticle permeability as well as the decrease in water loss and improvement in drought tolerance in wallflowers. GC‒MS analysis revealed that the wax components included alkanes, alcohols, aldehydes, and fatty acids. The most abundant wax compound in this plant was alkanes (50%), the most predominant of which was C29. The relative abundance of these compounds increased significantly in tetraploid plants under drought stress conditions. These findings revealed that tetraploid-induced wallflowers presented upregulation of multiple drought-related and wax biosynthesis genes; therefore, polyploidization has proved useful for improving plant drought tolerance.

Circadian rhythm response and its effect on photosynthetic characteristics of the Lhcb family genes in tea plant.

BACKGROUND: The circadian clock, also known as the circadian rhythm, is responsible for predicting daily and seasonal changes in the environment, and adjusting various physiological and developmental processes to the appropriate times during plant growth and development. The circadian clock controls the expression of the Lhcb gene, which encodes the chlorophyll a/b binding protein. However, the roles of the Lhcb gene in tea plant remain unclear. RESULTS: In this study, a total of 16 CsLhcb genes were identified based on the tea plant genome, which were distributed on 8 chromosomes of the tea plant. The promoter regions of CsLhcb genes have a variety of cis-acting elements including hormonal, abiotic stress responses and light response elements. The CsLhcb family genes are involved in the light response process in tea plant. The photosynthetic parameter of tea leaves showed rhythmic changes during the two photoperiod periods (48 h). Stomata are basically open during the day and closed at night. Real-time quantitative PCR results showed that most of the CsLhcb family genes were highly expressed during the day, but were less expressed at night. CONCLUSIONS: Results indicated that CsLhcb genes were involved in the circadian clock process of tea plant, it also provided potential references for further understanding of the function of CsLhcb gene family in tea plant.

Differential expression and analysis of extrachromosomal circular DNAs as serum biomarkers in pulmonary arterial hypertension.

BACKGROUND: Extrachromosomal circular DNAs (eccDNAs) have been reported to play a key role in the occurrence and development of various diseases. However, the characterization and role of eccDNAs in pulmonary arterial hypertension (PAH) remain unclear. METHODS: In the discovery cohort, we first explored eccDNA expression profiles by Circle-sequencing analysis. The candidate eccDNAs were validated by routine polymerase chain reaction (PCR), TOPO-TA cloning and Sanger sequencing. In the validation cohort, 30 patients with PAH and 10 healthy controls were recruited for qPCR amplification to detect the candidate eccDNAs. Datas at the baseline were collected, including clinical background, biochemical variables, echocardiography and hemodynamic factors. Receiver operating characteristic curve was used to investigate the diagnostic effect of the eccDNA. RESULTS: We identified a total of 21,741 eccDNAs in plasma samples of 3 IPAH patients and 3 individuals in good health, and the expression frequency, GC content, length distribution, and genome distribution of the eccDNAs were thoroughly characterized and analyzed. In the validation cohort, 687 eccDNAs were differentially expressed in patients with IPAH compared with healthy controls (screening threshold: |FC|≥2 and P < 0.05). Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed that the specific eccDNAs in IPAH were significantly enriched in calcium channel activity, the mitogen-activated protein kinase pathway, and the wnt signaling pathway. Verification queue found that the expression of eccDNA-chr2:131208878-131,424,362 in PAH was considerably higher than that in healthy controls and exhibited a high level of accuracy in predicting PAH with a sensitivity of 86.67% and a specificity of 90%. Furthermore, correlation analysis disclosed a significant association between serum eccDNA-chr2:131208878-131,424,362 and mean pulmonary artery pressure (mPAP) (r = 0.396, P = 0.03), 6 min walking distance (6MWD) (r = -0.399, P = 0.029), N-terminal pro-B-type natriuretic peptide (NT-proBNP) (r = 0.685, P < 0.001) and cardiac index (CI) (r = - 0.419, P = 0.021). CONCLUSIONS: This is the first study to identify and characterize eccDNAs in patients with PAH. We revealed that serum eccDNA-chr2:131208878-131,424,362 is significantly overexpressed and can be used in the diagnosis of PAH, indicating its potential as a novel non-invasive biomarker.